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BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates in general to camera calibration and more 
particularly to a method and a system for self-calibrating a wide field-of-view 
camera (such as a catadioptric camera) using a sequence of omni-directional 
images of a scene obtained from the camera. 

2. Related Art 

Cameras having a wide field of view are important in several diverse 
applications. By way of example, wide field-of-view cameras are used in 
applications such as surveillance, video conferencing, multimedia-based 
applications (such as games), autonomous navigation, advertising (especially in 
the areas of real estate and tourism) and scene reconstruction. Wide field-of- 
view cameras are used to capture most or all of a region of interest (or scene) 
while using a minimum number of cameras. 

One type of wide field-of-view camera is a catadioptric camera. In 
general, a catadioptric camera uses a combination of reflecting surfaces (such as 
mirrors) and lenses to achieve an omni-directional view of a scene. These 
reflecting surfaces, however, tend to cause image distortion. Image distortion is 
inherent in a catadioptric camera because light rays at a periphery of the field of 
view are severely curved to allow the use of a smaller, inexpensive image 
detectors (such as a charge-coupled device or CCD). Otherwise, prohibitively 
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large and expensive image detectors would be needed to capture the entire field of 
view. The image distortions tend to pull points toward an optical center of the 
camera, thus making straight lines within a scene appear curved in an image of the 
scene. 

5 If optical characteristics of a catadioptric camera system are known, then the 

parameters that map the curves in the image into straight lines can be applied. 
Unfortunately, however, these parameters are seldom known accurately and optical 
characteristics tend to vary between particular camera systems. This requires that a 
catadioptric camera be calibrated in order to determine calibration parameters that 
10 will characterize incident light rays and thus account for optical distortion. 

y3 Several types of calibration techniques are used to calibrate catadioptric 

in 

cameras. In one calibration technique the camera photographs a calibration pattern 
and uses the pattern to calibrate the camera. One disadvantage, however, of this 
technique is that the calibration pattern is large and using the pattern is awkward 
1 5 and time-consuming. Another calibration technique involves a user selecting a 
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number of points along straight lines of the scene (which are curved in the image). 
One disadvantage with this technique, however, is that the user must take the time 
to supply the calibration system with the data needed to correctly calibrate the 
0 camera. In reality, both of these calibration techniques may be difficult to implement 
20 in consumer-based products because it is unlikely a consumer would want to or 
could correctly perform such complicated and time-consuming calibration 
techniques. Moreover, if the camera calibration is performed incorrectly any results 
obtained from the camera will tend to be erroneous. 

Accordingly, there exists a need for a calibration method and system for a 
25 catadioptric camera that is accurate and simple to use. This calibration method and 
system should be capable of self-calibrating using images of a scene obtained from 
camera, thereby eliminating the need for calibration patterns or user-supplied points. 
In addition, this self-calibration method and system would provide a wide variety of 
calibration parameters to enable accurate results from the catadioptric camera. 
30 Whatever the merits of the above-mentioned systems and methods, they do not 
achieve the benefits of the present invention. 



2 



Docket No. 149508.1 

SUMMARY OF THE INVENTION 

To overcome the limitations in the prior art as described above and other 
limitations that will become apparent upon reading and understanding the present 
specification, the present invention includes a method and system for self-calibrating 
a wide field-of-view camera (such as a catadioptric camera) using a sequence of 
omni-directional images obtained from the camera. Based on characteristics of 
catadioptric imaging, the present invention self-calibrates the camera by determining 
a consistency of tracked features (such as points, lines and planes) across the 
sequence of images. This self-calibration method provides useful calibration 
parameters, including an aspect ratio and an image skew. 

The system and method of the present invention is accurate and simple to 
use. In particular, the present invention is self-calibrating and uses images from a 
particular scene environment to calibrate the camera, thereby eliminating during 
calibration the need for calibration patterns, user-selected points, knowledge of 
camera motion or knowledge of the scene geometry. Moreover, the present 
invention is capable of determining a number of calibration parameters that ensure 
an accurate and thorough calibration so that the catadioptric camera system can 
produce accurate images of a particular scene. 

In general, the calibration method of the present invention self-calibrates a 
catadioptric camera by obtaining a sequence of omni-directional images 
representing a scene, tracking features across the sequence and estimating a set of 
calibration parameters based on catadioptric imaging characteristics. The 
estimation of calibration parameters is performed by establishing a relationship 
(such as an objective function) between the feature tracks and the calibration 
parameters and estimating the calibration parameters using an error metric. An 
optimization technique is used to minimize the error metric and recover an optimal 
set of calibration parameters. The present invention uses at least two techniques to 
recover the calibration parameters, including an algebraic error metric and a 
preferred and more robust image error metric. In addition, the present invention 
generally determines more calibration parameters than other calibration techniques 
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including, for example, an aspect ratio and an image skew. The present invention 
also includes a system for self-calibrating a wide field-of-view camera using a 
sequence of omni-directional images that incorporates the method of the present 
invention. Moreover, the present invention includes a technique for reformulating a 
5 projection equation (such as used in the working example) such that the projection 
equation is equivalent to that of a typical rectilinear perspective camera. This allows 
analyses (such as structure from motion) to be applied in the same direct manner as 
for rectilinear image sequences. 

Other aspects and advantages of the present invention as well as a more 
10 complete understanding thereof will become apparent from the following detailed 

description, taken in conjunction with the accompanying drawings, illustrating by way 

Q 

y3 of example the principles of the invention. Moreover, it is intended that the scope of 

in 

t g the invention be limited by the claims and not by the preceding summary or the 
H following detailed description. 

m 15 

f BRIEF DESCRIPTION OF THE DRAWINGS 

£3 The present invention can be further understood by reference to the following 

Co 

|W description and attached drawings that illustrate the preferred embodiments. Other . 
jpS features and advantages will be apparent from the following detailed description of 
0 20 the invention, taken in conjunction with the accompanying drawings, which illustrate, 
by way of example, the principles of the present invention. 

Referring now to the drawings in which like reference numbers represent 
corresponding parts throughout: 

FIG. 1 is a block diagram illustrating an apparatus for carrying out the present 
25 invention. 

FIG. 2 is an overall block diagram of a catadioptric camera system 
incorporating the present invention. 

FIG. 3 is a general block diagram illustrating a method and system for self- 
calibrating a catadioptric camera using a sequence of omni-directional images in 
30 accordance with the present invention. 
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FIG. 4 is a general flow diagram of a self-calibration method of the present 
invention. 

FIG. 5 is a detailed flow diagram illustrating a working example of FIG. 4 of 
the present invention. 

FIG. 6 illustrates a cross-section of a portion of the catadioptric camera used 
in the working example; 

FIG. 7 illustrates the epipolar plane passing through a 3D point in space, two 
corresponding point features imaged at two catadioptric camera locations, as well as 
the two camera centers. 



DETAILED DESCRIPTION OF THE INVENTION 

In the following description of the invention, reference is made to the 
accompanying drawings, which form a part thereof, and in which is shown by way of 
illustration a specific example whereby the invention may be practiced. It is to be 
^ 15 understood that other embodiments may be utilized and structural changes may be 
made without departing from the scope of the present invention. 



I. Exemplary Operating Environment 



fU FIG. 1 and the following discussion are intended to provide a brief, general 

£3 

P 20 description of a suitable computing environment in which the invention may be 

implemented. Although not required, the invention will be described in the general 
context of computer-executable instructions (such as program modules) being 
executed by a computer. Generally, program modules include routines, programs, 
objects, components, data structures, etc. that perform particular tasks or implement 

25 particular abstract data types. Moreover, those skilled in the art will appreciate that 
the invention may be practiced with a variety of computer system configurations, 
including personal computers, server computers, hand-held devices, multiprocessor 
systems, microprocessor-based or programmable consumer electronics, network 
PCs, minicomputers, mainframe computers, and the like. The invention may also be 

30 practiced in distributed computing environments where tasks are performed by 

remote processing devices that are linked through a communications network. In a 
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distributed computing environment, program modules may be located on both local 
and remote computer storage media including memory storage devices. 

With reference to FIG. 1, an exemplary system for implementing the invention 
includes a general-purpose computing device in the form of a conventional personal 
5 computer 100, including a processing unit 102, a system memory 104, and a system 
bus 106 that couples various system components including the system memory 104 
to the processing unit 102. The system bus 106 may be any of several types of bus 
structures including a memory bus or memory controller, a peripheral bus, and a 
local bus using any of a variety of bus architectures. The system memory includes 
1 0 read only memory (ROM) 1 1 0 and random access memory (RAM) 1 1 2. A basic 
input/output system (BIOS) 1 14, containing the basic routines that help to transfer 

if i 

information between elements within the personal computer 100, such as during 
■Jjj start-up, is stored in ROM 110. The personal computer 100 further includes a hard 

H disk drive 1 16 for reading from and writing to a hard disk (not shown), a magnetic 

fftj 

1 5 disk drive 1 18 for reading from or writing to a removable magnetic disk 120, and an 

^ optical disk drive 122 for reading from or writing to a removable optical disk 124 

# 

O (such as a CD-ROM or other optical media). The hard disk drive 116, magnetic disk 

jjla, drive 128 and optical disk drive 122 are connected to the system bus 106 by a hard 

P J disk drive interface 126, a magnetic disk drive interface 128 and an optical disk drive 

£3 

P 20 interface 130, respectively. The drives and their associated computer-readable 

media provide nonvolatile storage of computer readable instructions, data structures, 
program modules and other data for the personal computer 100. 

Although the exemplary environment described herein employs a hard disk, a 
removable magnetic disk 120 and a removable optical disk 124, it should be 
25 appreciated by those skilled in the art that other types of computer readable media 
that can store data that is accessible by a computer, such as magnetic cassettes, 
flash memory cards, digital video disks, Bernoulli cartridges, random access 
memories (RAMs), read-only memories (ROMs), and the like, may also be used in 
the exemplary operating environment. 
30 A number of program modules may be stored on the hard disk, magnetic disk 

120, optical disk 124, ROM 110 or RAM 112, including an operating system 132, 
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one or more application programs 134, other program modules 136 and program 
data 138. A user (not shown) may enter commands and information into the 
personal computer 100 through input devices such as a keyboard 140 and a pointing 
device 142. In addition, a camera 143 (or other types of imaging devices) may be 
5 connected to the personal computer 100 as well as other input devices (not shown) 
including, for example, a microphone, joystick, game pad, satellite dish, scanner, or 
the like. These other input devices are often connected to the processing unit 102 
through a serial port interface 144 that is coupled to the system bus 106, but may be 
connected by other interfaces, such as a parallel port, a game port or a universal 
1 0 serial bus (USB). A monitor 146 or other type of display device is also connected to 
the system bus 106 via an interface, such as a video adapter 148. In addition to the 

D 

^ monitor 146, personal computers typically include other peripheral output devices 
Ml (not shown), such as speakers and printers. 

).£ The personal computer 100 may operate in a networked environment using 

q 15 logical connections to one or more remote computers, such as a remote computer 
H 150. The remote computer 150 may be another personal computer, a server, a 
p router, a network PC, a peer device or other common network node, and typically 
includes many or all of the elements described above relative to the personal 



fU computer 100, although only a memory storage device 152 has been illustrated in 
Q 



20 FIG. 1. The logical connections depicted in FIG. 1 include a local area network 

(LAN) 154 and a wide area network (WAN) 156. Such networking environments are 
commonplace in offices, enterprise-wide computer networks, intranets and the 
Internet. 

When used in a LAN networking environment, the personal computer 100 is 
25 connected to the local network 1 54 through a network interface or adapter 1 58. 

When used in a WAN networking environment, the personal computer 100 typically 
includes a modem 160 or other means for establishing communications over the 
wide area network 156, such as the Internet. The modem 160, which may be 
internal or external, is connected to the system bus 106 via the serial port interface 
30 144. In a networked environment, program modules depicted relative to the 

personal computer 100, or portions thereof, may be stored in the remote memory 
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storage device 152. It will be appreciated that the network connections shown are 
exemplary and other means of establishing a communications link between the 
computers may be used. 

5 II. Introduction 

The method and system of the present invention include self-calibrating a 
catadioptric camera system to obtain a set of calibration parameters. These 
calibration parameters are used, for example, to remove distortion in images 
captured using the catadioptric camera system. The present invention uses a 
1 0 sequence of images taken by the camera system of a scene (or environment) and 

„ uses this sequence of images to calibrate the camera system. The present 

M 

invention does not use special calibration patterns or need knowledge of camera 
t 2 motion or scene geometry. Instead, the present invention uses the consistency of 
^ features tracked across at least of portion of the image sequence and uses these 
£0 15 tracked features to find calibration parameters based on the characteristics of 

catadioptric imaging. The self-calibration method and system of the present 
p invention can obtain calibration parameters including, for example, aspect ratio and 

an image skew. The present invention also includes a method of transforming a 
p projection equation for a catadioptric camera system into an equivalent rectilinear 
£3 20 perspective using a transformation. 

III. General Overview 

As shown in FIGS. 2-7 for the purposes of illustration, the invention is 
embodied in a method and a system for self-calibrating a catadioptric camera 

25 system using a sequence of omni-directional images obtained from the camera 
system. FIG. 2 is an overall block diagram of a catadioptric camera system 200 
incorporating the present invention. The camera system 200 illustrated is one 
example of several camera systems that may incorporate the self-calibration method 
and system of the present invention. In general, the catadioptric camera system 200 

30 includes a catadioptric camera 205 including an image detector 210 (such as a 
charge-coupled device or CCD), a lens 215 (such as a telecentric lens) and a 
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reflecting surface 220 (such as a mirror) that is directed toward a region of interest or 
scene (not shown). Typically, the catadioptric camera 205 (including the image 
detector 210, lens 215, and reflecting surface 220) are located on a camera platform 
230 that permits the catadioptric camera 205 to be translated or rotated about at 
5 least some of the three-dimensional axes (not shown). A processor is connected to 
the catadioptric camera 205 and performs operations on data captured by the 
catadioptric camera 205 and calibration of the catadioptric camera system 200. 

Generally, self-calibration of the catadioptric camera system 200 occurs in 
accordance with the present invention by positioning the catadioptric camera 205 on 
10 a scene. The camera platform 230 may be rotated and translated as needed to 
position the catadioptric camera 205 such that the desired scene is captured. The 
scene does not need to contain any type of calibration chart (and the like) and a user 
^ does not need to select points within the scene for calibration to occur. Rather, the 
present invention self-calibrates the catadioptric camera system 200 by obtaining a 
1 5 sequence (or collection) of omni-directional images 250 of the scene. This 

sequence of omni-directional images 250 has not been processed by the processor 



O 240 and generally contains image distortions. The processor 240 receives the 

m 



sequence of omni-directional images 250 and sends them to a self-calibration 
module (not shown) to be used in the calibration process, which is discussed further 



20 below. 

After the catadioptric camera system 200 has been calibrated and calibration 
parameters obtained from the self-calibration module, the sequence of omni- 
directional images 250 obtained from the image detector 21 0 is sent to the 
processor. The processor applies the calibration parameters obtained from the self- 

25 calibration process of the present invention and outputs a collection of processed 
images 260. These processed images 260 may be used by other application 
modules, such as, for example, scene reconstruction, video conferencing or 
surveillance applications. 

FIG. 3 is a general block diagram illustrating a method and system for self- 

30 calibrating the catadioptric camera system 200 using the sequence of omni- 
directional images 250 in accordance with the present invention. In particular, the 
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image detector 210 obtains the sequence of omni-directional images 250 of the 
scene. This sequence of omni-directional images 250 includes several individual 
images 300 ranging from image(1) to image(N). The sequence of omni-directional 
images 250 is sent to a tracker 310 that detects features (such as points, lines and 
5 planes) within the sequence of omni-directional images 250 and tracks the positions 
of these features across the sequence 250. A self-calibration module 320 of the 
present invention receives the tracked features provided by the tracker 310 and, as 
explained in detail below, determines calibration parameters for the catadioptric 
camera system 200. These calibration parameters are sent as output for use by the 
1 0 catadioptric camera system 200 for providing accurate imaging results. 

u 

0 IV. General Operation of the Invention . 

FIG. 4 is a general flow diagram of a self-calibration method of the present 
H* invention. Generally, the self-calibration method determines calibration parameters 
P 15 by using a consistency of tracked features across a sequence of images based on 
H the characteristics of catadioptric imaging. These calibration parameters are used 
□ by the catadioptric camera system 200, for example, to eliminate distortion in images 
j|\ caused by the bending of light rays (i.e., to help make straight lines appear straight). 

fU The calibration process of the present invention is autonomous, and does not 
p 20 require user intervention such as special calibration patterns or selection of points on 
the image. Further, the present invention is convenient because no knowledge is 
needed of camera, motion or scene geometry. 

In particular, the self-calibration process of the present invention inputs the 
sequence of omni-directional images 250 (such as from the image detector 210). 
25 Features are tracked across the sequence of omni-directional images (box 410). 
Given an initial set of calibration parameters, the epipolar geometry, which links 
feature locations across two images within the sequence 250, can be defined (box 
420). Once the epipolar geometry has been characterized, an objective function can 
be defined as the deviation from the epipolar geometry for features across two 
30 images (box 430). In other words, the objective function is the error between the 
actual and predicted feature locations based on epipolar geometry. 
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Once the objective function is defined, optimal calibration parameters can 
then be extracted by employing an optimization technique that will optimize the 
objective function (box 440). This optimization may include, for example, in a 
preferred embodiment, minimizing the objective function. The output of the self- 
5 calibration module 320 is the result of using the optimization technique to optimize 
the objective function and is the optimal calibration parameters (box 450). These 
optimal calibration parameters can then be applied to a variety of applications such 
as 3D structure recovery using the catadioptric camera system 200. 

10 V. Details of the Invention and Working Example 

FIG. 5 is a detailed flow diagram illustrating a working example of FIG. 4 of 
D the present invention. This working example is provide for illustrative purposes and 
Iff is only one of several ways in which the present invention may be implemented. 
^ The structural and mathematical details of this working example are discussed in 

^15 detail below. 

CO 

^ Referring to FIG. 5, a sequence of omni-directional images obtained from the 

catadioptric camera system is received as input (box 500). Next, point features are 
tracked across this sequence of images (box 510). Correspondences between point 
features across two different images are identified and used in the objective function 
It 20 to recover the calibration parameters (box 520). An objective function is defined as 
an error metric defined in terms of the deviation of pairwise correspondences from 
the epipolar geometry (box 530). The present invention includes two error metrics, 
both of which are defined in detail below. A preferred error metric is an image-based 
error metric that defines an image distance to an epipolar curve. The image-based 
25 error metric is preferred because it is more robust. An alternative error metric is an 
algebraic error metric that is not as robust as the preferred image-based error 
metric. 

' The error metric is used to find an optimal set of calibration parameters that 
minimizes the error of the objective function. The optimum set of calibration 
30 parameters can be obtained by minimizing the objective function using an 

optimization technique (box 540). These optimization techniques are known to 
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those skilled in the art. Prior to using the optimization technique, initial values of at 
least some of the unknown calibration parameters are selected as a starting point for 
the optimization (box 550). These initial values may be reasonable guesses or 
based on another calibration technique such as, for example, a circle-based 
calibration technique discussed in detail below. Optimal calibration parameters then 
are obtained from this optimization (box 560). 

In this working example, the catadioptric camera system is a stand alone, 
movable system that is capable of capturing long sequences of omni-directional 
images. One application of this system is to reconstruct complex large 
environments, such as an entire floor of a building, from the captured images only. 
The camera system includes a color camera having a resolution of 1 140 *1030 
pixels and a catadioptric attachment including a panoramic lens and paraboloid 
mirror. The camera and catadioptric attachment are located on a platform that 
includes a direct current (DC) motor that aids in positioning the camera. An 
uninterruptable power supply (UPS) is used to provide remote power, and 
processing is performed by a personal computer having a 450 MHz processor with 2 
gigabytes (GB) of random access memory (RAM) and a flat-panel display. 

The large RAM size enables slightly more than 1,500 omni-directional images 
to be captured on the fly in the raw image mode of 1 byte per pixel. Specifically, up 
to 1,500 images at 6.7 Hertz (Hz) and a resolution of 1 1 140 x 1030 may be captured 
by this system. The rate of capture can be as high as 1 1 .75 frames per second. 
With the use of the UPS, the system can perform standalone capture for about half 
an hour. 

In this working example, collections of images were taken under two 
conditions. The first condition was while the camera was being rotated about a 
vertical axis by the motor (which has variable speeds) while the platform was 
stationary. The other condition was while the camera was stationary with respect to 
the platform and the platform was manually moved. 

In the self-calibration method of the present invention embodied in this 
working example, point features were tracked across a sequence of omni-directional 
images. The camera system contained a set of "blind spots" caused by obstruction 
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of view by the supports for the mirror and camera. Point features tracked along the 
perimeter of these "blind spots" could not be assumed to be stationary across the 
collection, because it was possible that erroneous point feature tracks could occur 
due to occluded edges giving rise to spurious corners. In order to avoid the 
problems caused by the "blind spots " within the collection of images, a mask was 
created that was used by the self-calibration method to ignore points inside or within 
five pixels of the mask boundary. 

The point feature tracks in this working example were generated using a 
tracker developed by J. Shi and C. Tomasi entitled 'Good Features to Track" in 
CVPR, pages 593-600, Seattle, Washington (June 1994), the entire contents of 
which are hereby incorporated by reference. In this working example, the point 
features were not tracked across a long collection of images because of the highly 
distorting imaging characteristics of the catadioptric camera. Instead, a point feature 
tracking was performed across a collection of 20-30 images spanning between 15 
and 25 degrees. 

Imaging Characteristics of a Catadioptric Camera System 

In the self-calibration method of the present invention, a set of unknown 
calibration parameters are defined using imaging characteristics of the catadioptric 
camera system. The mathematical details of the present invention are presented 
below. 

The catadioptric camera system consists of a telecentric lens that is designed 
to be orthographic and a paraboloid mirror. FIG. 6 illustrates a cross-section of the 
image plane-mirror representation through the center (or focus) of the paraboloid 
mirror. 

The expression for the mirror surface is 



An actual observed image pixel location (w, v) is linked to a camera pixel location 
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(h',v') by a calibration matrix M: 
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where a is an aspect ratio, s is an image skew, and {p x ,p y ) is a principal point. 

The present invention also includes a method of transforming a projection 
equation for a catadioptric camera system into an equivalent rectilinear perspective 
using a transformation. In other words, the projection equation of the catadioptric 
camera system can be converted, using the following method, to a form of a 
rectilinear perspective camera through a transformation of the z coordinate. In 

particular, a point ? = (x,y,zf relative to O in FIG. 6 is mapped to camera pixel 

{u'yf . The actual 3D point on the mirror surface that P is projected onto is given 
by the projection equation 
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v' 



h 2 -{u' 2 +v' 2 ) 
2h 



(3) 



with A > 0 . Hence 



A 2 -(AV+AV) 
2h 



(4) 



Solving for X , gives 
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X = h 



^x 2 + y 2 + z 2 -z 
- - 



x 2 +y 2 +z 2 + z 



x 2 +y 2 +z 2 -hz 



(5) 



4x 2 + y 2 + z 2 +z |P| + z z' 



As a result, the projection equation can be rewritten as 
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Substituting Equation 7 into Equation 2, gives 
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This is the form for rectilinear perspective projection. Therefore, using this 
formulation, structure from motion on the tracks in the omni-directional image 
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sequence can be applied in the exact manner as for rectilinear image sequences. 
The actual 3D points (x,y,zf can then be computed from the resulting "pseudo" 3D 
points (x,y,z') T using the relationship 



since z' = ^]x 2 +y 2 + z 2 + z . This projection equation transformation of the present 
invention can be used in the determination of calibration parameters to recover 
structure from image collections. 

Epipolar geometry may also be considered when determining the calibration 
parameters. In particular, for a pair of images with indices 1 and 2, corresponding 
points (in homogeneous coordinates to indicate their positions on the mirror surface) 
satisfy the epipolar constraint 



where E is the essential matrix associated with the image pair, and q jk is defined in 
Equation 3. 

FIG. 7 illustrates an epipolar plane for cameras positioned at Oi and 0 2 , 
wherein % and q 2 are the respective projections of 3D point Pj. From FIG. 7, the 
normal to the plane passing through both camera centers (with respect to 0 2 ) is 




qj 2 £q yl =0, 



(10) 



CD 



Hence, the equation of the plane is of the form 
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n 2 -p = 0oraw' + &v' + cw/ = 0 (l2) 
After substituting Equation 1 into Equation 12 and rearranging terms, 

u >2 +v >2- 2 ^u'-2^v'-h 2 =0 (13) 



This is the equation of a circle centered at 
Equation 13 may be rewritten as 
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From Equation 2, the epipOlar curve is obtained in the second image corresponding 
to q ;1 in the first image as 
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Circle-Based Calibration Technique 

The circle-based technique is used in this working example to determine 
initial values for some of the unknown calibration parameters. This circle-based 
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technique uses a single image to identify a bounding circle of the omni-directional 
image. This may be done manually or automatically by using a predefined 
threshold, finding a boundary and fitting a circle to the resulting boundary. The 
center of the circle is defined as an image principal point. Since the field of view of 
the omni-directional image is known (105° in the catadioptric attachment of this 
working example), a camera parabolic parameter can then be directly computed 
using the radius of the circle. 

Referring to FIG. 6, the relationship between the mirror parameter H and the 
radius of the omni-directional image may be determined. In particular, from 

h 2 —r 2 r 

z = andr = tan# = - (l6) 

2h z v ' 




gives 



'-it? ° 7) 



After some manipulation, the following expression is obtained 



h=r i±l±<L =rcot 0 (18) 

t 2 v ' 



Therefore, h may be computed if r = r max (i.e., the radius of the omni-directional 

image) is known corresponding to 6 = 0 max =1O5°. As explained further below, this 

circle-based calibration technique may be used to determine initial values for some 
of the unknown calibration parameters. 
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Estimation of Calibration Parameters 

The estimation of a set of unknown calibration parameters (i.e.,h 9 (p x9 p y ) , r , 

and s as defined in Equations 1 and 2) is accomplished by using a least-median 
error metric. In this working example, an exhaustive combination of pairs of images 
were used that were at least four frames apart and had a minimum of 10 
correspondences in an objective function to be minimized. This maximized the use 
of point features and avoided possible degenerations associated with very small 
camera motion. Estimation of essential matrix of pairs of images was performed 
using a method described by Z. Zhang entitled "Motion and structure from two 
perspective views: From essential parameters to Euclidean motion via fundamental 
matrix" in J. of the Optical Soc. of America A, 14(1 1):2938-2950 (1997), the entire 
contents of which are hereby incorporated by reference. 

To recover the calibration parameters, an objective function 



was defined where 7V pairs is the number of different image pairs, "med "refers to the 
median of the series of errors, E^, b i and e, are the image numbers corresponding 
to the zth pair, and S(p iy e) is the set of indices of point feature tracks that spans 
across at least images b i and e i . Ey can be an algebraic error metric 



pairs 



(19) 



tf = fe^J ! (20) 



or an image error metric 



EP = d^,E iq J + d^ A ,Ej q J (21) 



19 



Docket No. 149508.1 

where rf(m, n) is the image distance from m to the epipolar curve specified by n of 
the form specified in Equation 15. In this working example, the image-based error 
metric was used. 

In the present invention, estimation of the unknown calibration parameters 
5 (h, (p x ,p y ), a, and s) is performed using an optimization method. These 

optimization methods are known to those skilled in the art. In this working example, 
the Nelder-Mead simplex search algorithm was used. The initial values of h and 
(Px'Py) were determined by using a circle based calibration technique described 

above. The initial values of a and s were selected to as 1 and 0, respectively. This 
1 0 technique of estimating the unknown calibration parameters was used for this 
working example because the technique is guaranteed to converge to a local 
minimum and is simple to implement. 

Results 

15 The self-calibration method of the present invention as implemented in this 

working example were applied to a number of different image collections and 
produced results that were consistently better than existing calibration techniques. 
The known straight lines in a processed (or dewarped) images produced with the 
calibration parameters determined by the present invention appear to be straight. 

20 While existing calibration techniques did on occasion produce reasonable results, 
they were not as consistent as those of the self-calibration method using the 
preferred image-based error metric. Table 1 numerically illustrates that the self- 
calibration method of the present invention using the preferred image-based error 
metric produced the most favorable results. The numerical outputs shown in Table 1 

25 are for the same collection of images (both h and p are in pixels). 
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Circle-based 


Algebraic error 


Image error 


h 


345.33 


530.85 


408.28 


p 


(630.54, 


(609.36, 


(575.22, 




480.78) 


410.50) 


453.62) 


a 




0.988621 


0.952877 


s 




-0.001724 


0.002096 



Table 1 



C3 In order to verify that the calibration parameters obtained from the self- 

In calibration method using the preferred image-based error metric are superior, a 

it ^ validation was performed using calibration parameters obtained from the circle- 
Si based calibration technique, the self-calibration method using the algebraic error 

ESS 

jTT metric and the image-based error metric. Specifically, a mosaic of an office was 

* & extracted using a sequence of parallax-free images. The camera was rotated such 

o 

m that the rotation axis passed through the virtual projection center of the camera. 

1 0 This was accomplished by continually adjusting the position of the camera and 
O checking for parallax by rotating the camera and viewing the images. In this working 

example, a collection of 180 images were used, each about 1° apart. Because the 
maximum angular field of view for each camera was approximately 210°, the total 
span was about 389°, which was enough to wrap around itself. In addition, there 
1 5 was no enforcement of the closed loop constraint in constructing the mosaic. 

Using the calibration parameters obtained from the circle-based calibration 
technique, the mosaic appeared blurred. In contrast, the calibration parameters 
obtained using the other two techniques appeared more in focus. However, the 
calibration parameters obtained using the algebraic error metric did not close 
20 properly, probably due to underestimation of the rotation angle. The self-calibration 
method using the image-based error metric provided the best dewarped image. 
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The foregoing description of the preferred embodiments of the invention has 
been presented for the purposes of illustration and description. It is not intended to 
be exhaustive or to limit the invention to the precise form disclosed. Many 
modifications and variations are possible in light of the above teaching. It is 
intended that the scope of the invention be limited not by this detailed description of 
the invention, but rather by the claims appended hereto. 
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